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Abstract 

We propose a new method for exact analytical calculation of the ground-state energy of the 
Ising spin glass on strips. An outstanding advantage of this method over the numerical transfer 
matrix technique is that the energy is obtained for complex values of the probability describing 
quenched randomness. We study the ±J and the site-random models using this method for 
strips of various sizes up to 5 x oo. The ground-state energy of these models is found to have 
singular points in the complex-probability plane, reminiscent of Lee- Yang zeros in the complex- 
field plane for the Ising ferromagnet. The ±J Ising model has a series of singularities which 
may approach a limiting point around p ~ 0.9 on the real axis in the limit of infinite width. 

1 Introduction 

In spite of long history of intensive studies, there are few exact results known on models of 
spin glasses. To summarize analytical calculations of the ground-state properties of the Ising 
model, Derrida et al. [[IJ investigated linear chains, double-chain strips and the 3x3 square 
lattice systems using a variant of the transfer-matrix method to obtain exact expressions of 
the ground-state energy. Morita and Horiguchi || and Fechner and Blaszyk J§| investigated the 
linear triangular chain by summing up all the configurations of frustrated clusters. They derived 
analytical expressions of the ground-state energy for these systems. However, their methods 
are essentially an enumeration of possible configurations and therefore are very difficult to be 
generalized to systems with larger widths. 

Recently, interest in this problem has been renewed by Dress et al. who applied the method 
of Derrida et al. jij to the problem of error-correcting codes. In the present paper, we generalize 
their method to the calculation of the ground-state energy of the ± J model and the site-random 
model on strips of various width. An important reason to consider the latter site-random model 
is that almost all experiments of spin glasses are carried out on site-random systems. Although 
bond-random and site-random models would share many important properties, there should 
exist some distinctions, a part of which is clarified in the present paper. An advantage of this 
method over the enumeration methods is that ground-state energy of strips of various width with 
various kinds of randomness can be calculated systematically without elaborate enumeration of 
configurations. 

A useful aspect of the method used in the present paper is that the ground-state energy is 
obtained for general complex values of parameters to control probabilities of quenched random 
variables. This is in sharp contrast to the numerical transfer matrix technique which is applicable 
only to real probabilities. As long as the width of a strip is finite, the ground-state energy of 
models on the strip is analytic as a function of the probability of quenched randomness if 
this probability is restricted to real values. However, when the variable to control probability 
is analytically continued into the complex plane, the ground-state energy is found to have 
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singularities. This analytic continuation into the complex plane is analogous to the idea of Lee- 
Yang zeros for complex magnetic fields. Thus, a hint for the singularity of the ground-state 
energy of the two-dimensional spin glass may be obtained from the behavior of the singularities 
in the complex plane as the width of strip is increased. 

In §2, we explain the formulation of the present method using examples of the width- two 
(ladder) systems. Explicit formulas of the ground-state energy of the bond-random and site- 
random models are given. Investigation of systems with widths larger than two needs additional 
techniques, which are explained in §3. The exact ground-state energy is given for the case of 
width three. We could obtain locations of singularities of the ground-state energy only in 
limited regions of the complex probability plane when the width is four and five. This is due 
to limitations of computational resources (memory and CPU time). We should stress here that 
exact analytical expressions, not numerical values, of the ground-state energy can be obtained 
if sufficient computational power is available. In §4, behavior of singularities in the complex 
probability plane is discussed. Differences between bond-random and site-random models are 
described and their physical implications are discussed. 



2 Formulation of Exact Calculations 



In this section, we formulate the procedure for the evaluation of the exact ground-state 
energy of the ±J and the site-random models. The width-two (ladder) systems are treated as 
examples. This formulation is based on the method proposed by Derrida et al. Q and Dress et 
al. The main problem is how to evaluate averages over the distribution of randomness. 

2.1 ±J model 



by 



First, we consider the ± J Ising model on a ladder shown in Fig. [[]. The Hamiltonian is given 



H — — y^^KjOi-KTj + K' i a' i _ 1 a' i + Licncr'i) , 



(1) 



where the variable o-i(a' i ) = ±1 represents an Ising spin on the upper (lower) line, and the 
exchange interactions Ki, K[ and Li take the values +1 or —1 randomly with the probabilities 
p or 1 — p, respectively. We calculate the ground-state energy per spin in the limit n — » oo. 

As is well-known, transfer-matrix calculation of the partition function of Ising models is 
based on the constrained partition functions Z„ , in which the boundary spins are fixed as 
°"n,On = ±1 in the 2 x n-spin ladder. These constrained partition functions obey the following 
recursion relation 



(2) 
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using the transfer matrix T n given by 
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(3) 



where 



(4) 
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In the zero-temperature limit z S> 1, we have only to consider the leading term of z, namely, 



7+~ 
7- H 



A+- 

^4 77. 



y x n +2a n 
7 x n +2b n 



(5) 



The exponent x n represents the energy of the strip of length n with the boundary condition 
a — a' — +1. The exponents a n , b„ and c n correspond to the differences of energy between the 
systems with different fixed boundary conditions. 

Since we are interested only in the energy, we concentrate on the exponents x n , a n , b n and 
On. For the asymptotic form the recursion relation (^) is reduced to 



Xn + l 



= K + K' + L + x n + 2M, 

— max(—K' — L, —L + a n , —K — K — L + b n 
—K — L + Cn) — M, 

— m&x(—K — L, —K — K' — L + a n , 
—L + b n , —K' — L + c„) — M, 

= max(— K — K , — K -+- a n , 

-K' + b n ,Cn)-M, 



(6a) 
(6b) 
(6c) 
(6d) 

max(0, —K' + On, —K + b n , -K - K' + c n ) , (7) 

and K = isT n +i) K' = K' n+1 and L = L n +i- Because of the global reflection symmetry in 
the spin space, b n is equal to a n and c n is equal to zero. Then, the four recursion relations in 
(|6a|)-(|6d|) are reduced to the following two recursion relations 

x n+1 = K + K' + L + x n + 2M', (8a) 
a.n+1 — m&x(—K' — L,—L + a n ,—K — K' — L + a n , 



with 



M 



-K-L) — M', 



with 



M' = max(0, —K' + o n , -K + o„, —K - K') 
The ground-state energy per spin is expressed as 



(8b) 
(9) 
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n— >oo In V p 



lim — [x n + max(0, 2a„, 2a n , 0)] 

n — >oo ZTL 



lim — 

77 — >00 2n 



^ (im+i - im) + + max(0, 2a„) 



(10) 



where the double brackets <C ■ ■ ■ 3> represent the average over the bond configuration. The 
quantity — (x n +i — x n ) corresponds to the increase of the energy by the transfer from Z n to 
Z n +i- Equation ( paj ) means that the quantity (x n +i — x n ) is a function of a n , K n +i, K' n+1 and 
i n +i- Since K n +i, K' n+1 and L n +i are independent of each other and of a n , the probability 



P(a n ,K, 



Ti+i , K n+1 , Ln+i 



is decomposed as P(a„) ■ P(K n +i) ■ P(K n +i) ■ P(K n +i). Among the 



four factors in this product, only the distribution P{a n ) is nontrivial. The function P(a n ) can 
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be calculated if we analyze the change of state from the nth stage of transfer to the (n + l)th 
stage. Details are explained in Appendix The result is 

i?(p) = - 3 7 3p+23p r 40p 3 3+2 t 4 . (id 

y ' 2(1 + p + 3p 2 - 8p 3 + 4p 4 ) ^ ' 

This quantity is plotted as a function of real p between and 1 in Fig. ^. Note that this function 
is symmetric for the exchange of p with 1 — p as it should be. 

2.2 Site-random model 



Next, we consider the site-random Ising model on a ladder. The Hamiltonian is given by 
the same form as in eq. (Q), except for the expressions of the exchange interactions K n , K'„ 
and L n . The site-random model is composed of A and B ions with concentrations c and 1 — c, 
respectively. The value of the exchange interaction is +1 for the neighboring ion pair A- A and 
is -1 for A-B and B-B. 

Calculations of the ground-state energy in this model can be carried out quite similarly to 
those of the ± J model. Formulas ([!]) to ( Jioj ) do not need modifications. Only the evaluation of 
the final expression ([n^) requires different treatments as explained below. 

In the site-random model, the variables K n +i, K' n+1 and L n +i are functions of the type 
of ions X n and X' n which are either A or B in the nth row and ions X n +i and X n+1 in the 
(n + l)th row. Hence, the quantity (x n +i — x„) is a function of a„, X„, X' n , X n +i and X' n+1 . 
Let us remember that the variables X n +\ and X' n+1 located at the edge of the strip are random. 
That is, 

p m = {' ! « ! °«* = A b • <12) 

where X stands for X n +i or X' n+l . On the other hand, the variables X n and X' n are not random, 
since they can be regarded as having already been fixed in the previous step of transfer. Then, 
the distribution function is written as 

P(a n , X n , X n , X n +i, X n+1 ) 

= P(a n ,X n ,X' n )-P(X n+1 )-P(X' n+1 ) . (13) 

The non-trivial part of distribution function is P(a n , X n , X' n ). Evaluation of this P(a n , Xn, X' n ) 
proceeds quite similarly as in the case of the ± J model. Details are found in Appendix |Aj. Wc 
obtain the explicit expression of the ground-state energy as 

E(c) = - 3(1 - 3c 2 + 6c 3 -3c 4 ) 
K ' 2(l + C -c 2 )(l- C + c 2 ) • K > 

The value of this energy is drawn as a function of real c in Fig. ^. It is seen that this quantity 
is symmetric under the exchange of c and 1 — c as was discussed generally in Ref . Q . 

3 Systems with Larger Widths 



In the present section, we generalize the method explained in the previous section to strips 
with width lager than two for the ±J and the site-random models. 



4 



In the case of width 3, the constrained partition functions in the low temperature limit is 
expressed as 



(15) 



where the allocation of the exponents x n ,a n , b n and c n takes the global spin- flip symmetry into 
account. These exponents satisfy the following recursion relations analogous to eqs. ([Ba|) , (pb|) : 
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K + K' + K" + L + L" + x n + 2M", (16a) 
max(-K" - L 1 , -L' + a, —K' - K" - L' + b, -K' - L' +c, 
—K — K" — L + c, -K-L' + b, -K -K' + a, 

—K — K' — K") - M", (16b) 

m&x(-K' -L-L', —K' - K" - L - L' + a,-L - L' + b, 

—K" - L - L' + c,-K - K' - L - L' + c,-K - K' -K" -L-L' + b, 

-K -L-L' + a, —K — K" — L — L") - M", (16c) 

m&x(-K' - K" - L, -K' - L + a, -K" - L + b,-L + c, 
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: max(0, -K" + a, -K + b, —K' - K" + c,-K + c, -K - K" + b, 

—K — K' + a, —K — K' — K") , (17) 

- b n , C = C n , K — K n + 1, K' = K' n +1, K" — K" +1 , L — Ln+1 Slid L' = L'n + l- 

We apply the above method to the ±J and the site-random models. However, as was 
mentioned in Appendix [a], the same approach to the case of width 2 requires a symbolic inverse 
of a matrix of size 28 x 28 or 60 x 60, which requires exceeding amount of CPU time. The 
following trick makes it possible to calculate the exact ground-state energy without explicit 
symbolic inverse of large matrices. 

Equation (pSl|) can be written as 



with 



where a 



M" 



P(a„) = (S 1 )i,N = 



det(S) 



(18) 



where the index i (= 1, • ■ ■ , N) of Pi stands for an element of the set {a n } and Sn,i denotes the 
co-factor of S. Then, eq. ([To]) is expressed explicitly as 



E(p) 



-Z^2 Pi ■ (Zn+l 



2det(5) 



n+1 



3>n)i 



(19) 



Let us recall here that Sn,% and det(S') are polynomials of p with integer coefficients. Since 
it is easy to calculate Sn,i or det(S) for given p numerically, we can determine the coefficients 
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of the polynomials from numerical values at several points: If the degree of the polynomial is 
m, numerical values of the polynomial at (m + 1) points are sufficient to determine the exact 
values of coefficients. The same technique holds for the site-random model. 

In this way we have carried out the average of eq. ( |ici| ) and obtained the ground-state energy 

as 



with 



N(p) = 5 - 25p + 158/ - 596/ + 1703/ - 3774/ 
+6322/ - 8368/ + 9892/ - 11536/ 
+12304/° - 10112/ 1 + 5568/ 2 

-1792/ 3 + 256/ 4 , (21) 
D{p) = 3(1 - 3p + 22/ - 82/ + 231p 4 
-516/ + 844/ - 960/ 

+720p 8 - 320/ + 64/°) (22) 



in the ±J model, and 



with 



E(C) -~DMC)' (23) 



N S r(c) = 5 - 15c 2 + 22c 3 + 19c 4 - 52c 5 + 40c 6 
-112c 7 + 370c 8 - 580c 9 + 468c 10 

-192c 11 + 32c 12 , (24) 
D S r{c) = 3(1 - c + 3c 2 - 4c 3 + 6c 4 - 6c 5 + 2c 6 ) 

x (1 + c - 3c 2 + 6c 4 - 6c 5 + 2c 6 ) (25) 

in the site-random model. These results @ and (f23|) are plotted as functions of p or c in Figs, a 
and |E|. 

In the case of width 4, the number of the bases of the transition matrix S defined in Ap- 
pendix |a| exceeds two hundred. Accordingly the degrees of polynomials of the denominator and 
numerator of the expression of the energy ( |l9| ) are much larger than in the case of width 3. 
Thus we have to determine the coefficients of polynomials of very large degrees from numerical 
values of the polynomials at many points. This procedure requires prohibitively high- precision 
calculations, which we were not able to perform. 

Instead, we determined the locations of singularities (where det(S) = 0) in restricted regions 
of the complex probability plane. First, the value of det(S) is scanned in the complex plane 
with the interval of 0.01 for real and imaginary parts of p or c. Zeros of det(S') correspond 
to the crossing points of two sets of curves, the first one corresponding to Re[det(S')] = and 
the second to Im[det(S)] = 0. These curves are displayed in Fig. m. Thus, the regions in 
which zeros are likely to be located are specified from the crossing points. Next, the zeros in 
the specified regions are precisely evaluated numerically using the subroutine of finding zeros 
of a function in the complex plane. In some cases, zeros of det(S') coincide with zeros of the 
numerator of eq. ([l9|). In such cases the energy is not singular. This type of zeros of det(S) 
should be excluded from our considerations. Direct numerical evaluation of the energy at such 
points reveals whether or not a zero of det(S) is a real singularity of the energy. We have in 
this way singled out real singularities of the energy. Note that the present method gives much 
more precise results than directly locating divergences of the energy in the complex plane which 
accompanies large rounding errors. 
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In the present analyses, we have concentrated on the region —0.5 < Re(p, c) < 1.5 and 
— 1 < Im(p, c) < 1 due to the limit of computational time. Although some singularities lie outside 
this region, physically important ones are located near the real axis in the range < p, c < 1. 
Therefore, for the purpose of investigation of effects of singularities in real physical systems, it 
is sufficient to concentrate our attention to the above mentioned region in the complex plane. 

In the case of width 5, we could locate only one singularity nearest to the range < p < 1 
for the ±J model by using the same method as the case of width 4. The site-random model 
was not investigated for the width 5. 

Let us note that the unit of strips is not limited to squares. For example, a strip composed 
of triangles can be treated, as shown explicitly for the case of width 2 in Appendix IbI 



4 Summary and Discussion 



In the present paper, we have generalized the method for the calculation of the exact ground- 
state energy proposed by Derrida et al. |lj and Dress et al. Q|. The ground-state energy of 
finite-width strips can be obtained by using this method. 

We have applied the method to the ± J Ising model on strips with width 2 to 5 and the site- 
random Ising model with width 2 to 4. In the case of widths 2 and 3, the exact ground-state 
energy and the singularities in the complex-probability plane of randomness have been obtained 
explicitly. In the case of width 4 and 5, the exact location of singularities near the physical 
region (0 < p, c < 1) have been identified by numerical calculation. 

In the ±J model, Fig. ^(a), there seems to exist a series of singularities which may approach 
a limiting point around p ~ 0.9 on the real axis as the size of the strip grows. In the site-random 
model, Fig. ^(b), such series does not seem to exist within the present calculations. These results 
suggest that the ground-state energy of the ±J model may have a singularity at the critical 
concentration p c ~ 0.89, [Q] and that the ground-state energy of the site-random model does 
not have any singularity at the critical concentration c c ~ 0.63 jf|. Further investigations are 
necessary to clarify the significance of these observations. 
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A Evaluation of the Functions P(a n ) and P(a n , X n , X' n ) 



T he recursion relation (^) has a stationary distribution P(a n ). As has been mentioned in 



2.1 



the explicit form of this distribution function is necessary for exact calculation of the 
ground-state energy. It is found easily by explicit iteration of the recursion relation ( jjb| ) start- 
ing from various finite values of a n that the set {a n } is restricted to a closed finite one. That 
is, we find {a n } — {2, 1,0, —1, —2} in the ±J model of width 2. Similarly, in the site-random 
model with width 2, it is necessary to obtain the joint distribution function P(a n , X„, X' n ) 
as explained in § UJA. It is again found by explicit iteration that the set {a n , X n , X' n } con- 
sists of (0,A,A), (-1,A,A), (-2,A,A), (0,A,B), (— 1, A, B), (-2,A,B), (0,B,A), (-1,-B^A), 
( — 2,B,A), (2, _B, B), (1, B, B) and (0, B,B) as the stationary state of the recursion relation (pb|). 
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Bearing these observations in mind, we construct a matrix S which transforms a n to a n +i- 
In the case of the ± J model with width 2, we find from eq. (|8q) 



\a„+i = +2) |on+i = +1) \a n +i = 0) \a n +i = —1) |ot»+i = —2} 



S = 



p 2 q + q 3 p 2 q + q 3 





p 2 q + 2pq 2 + q 3 



2pq 2 



2pq 2 

p 2 q + 2pq^ + g a 

p 3 + 3pq 2 p 3 + 3pq 2 3p 2 q + q 3 3p 2 q + q 3 

p 3 + 2p 2 q + pq 2 



2p 2 q 



2p 2 q 



p 3 + pq 2 p 3 + pq 2 



| On = + 2 ) 
|«n = +1) 

|o« = 0) , (26) 
|a« = -1) 
a n = -2) 



with g = 1 — p. The stationary distribution P(a n ) is given by the eigenvector of S with the 
eigenvalue unity: 

/ P{2) \ ( P(2) 
= 8 



P(-2) 



where P(a n ) satisfies 



J] p K) = 1 

{<*„} 



It is convenient to introduce the matrix S defined by 



(27) 
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where N denotes the size of the set {a n } (N = 
(pgj) are written in a single relation 


5 in 


the case of eq. (pj|)). Then eqs. 
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(30) 



Thus, we obtain 
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(31) 



Similarly, in the case of the site-random model with width 2, the distribution function can 
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be obtained by usinj 



g the transition matrix 
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The matrix S or S depends on models. The size of the matrix increases exponentially as 
the system size increases. The sizes of the matrices are 5, 28, 286, 3400 for widths 2, 3, 4, 5 in 
the ±J model, and 12, 60, 528 for widths 2, 3, 4 in the site-random model. Thus the explicit 
evaluation of an inverse matrix by symbolic manipulation becomes difficult. 



B Strip Consisting of Triangles 



We have calculated the ground-state energy of the one- dimensional ±J model on the linear 
triangular chain (Fig. pi) by using the same method. The ground-state energy is given by 



2(2 - 2p + 6p 2 - 4p 3 ^ 
1 + 3p 2 + 2p 3 

This model has already been "solved" by Fechner and Blaszyk Q . Although this expression 
is consistent with their Fig. 3, their expression of the ground-state energy is neither consistent 
with this expression nor with their own figure. 
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Figure Captions 



Figure 1: The ladder of the ± J model with random ferromagnetic (solid line) and antiferromagnetic 
(dashed line) interactions. 

Figure 2: The ground-state energy of the ladder system of the ±J model as a function of real p. 

Figure 3: The ground-state energy of the ladder system of the site-random model as a function of 
real c. 

Figure 4: The ground-state energy of the ± J model on the 3 x oo-spin strip. 

Figure 5: The ground-state energy of the site-random model on the 3 x oo-spin strip. 

Figure 6: Curves on which the real part of det(S') vanishes, (a) and (c), and the imaginary part 
vanishes, (b) and (d), for the 4 x co-spin strips. The figures (a) and (b) are for the ±J model, and 
(c) and (d) are for the site-random model. Only the quarter of the whole region (Re(p, c) > 0.5, 
Im(p, c) > 0) is shown, because the models are symmetric under the exchange p, c «-> p* , c* (complex 
conjugate) or p, c <-> 1 — p, 1 — c. 

Figure 7: Singularities of the ground-state energy in the ± J model (a) for the 2 x oo- to 5 x co-spin 
strips, and in the site-random model (b) for the 2 x oo- to 4 x oo-spin strips. The symbols o, □, 
O and A correspond to the singularities of the 2 x oo-, 3 x oo-, 4 x oo-, and 5 x oo-spin strips, 
respectively. The symbol • in (a) denotes the transition point from the ferromagnetic phase to a 
non-ferromagnetic phase along the real axis for the ±J model, p c = 0.89 ±0.01 0. Only the quarter 
of the whole region is shown because of the symmetry as explained in relation to Fig. ^|. 

Figure 8: The linear triangular chain with random ferromagnetic (solid line) and antiferromagnetic 
(dashed line) interactions. 
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